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We show that the new measurement of the effective number of neutrinos (iV c ff) by the Planck 
experiment can be used to set a robust lower bound on the mass of cold thermal dark matter 
of O(MeV). Our limit applies if the dark matter remains in thermal equilibrium by coupling to 
electrons and photons or through interactions with neutrinos, and applies regardless of whether the 
dark matter annihilation cross-section is s-wave or p-wave. To illustrate our bounds we apply them 
to a model of a supersymmetric neutralino annihilating to neutrinos, via a light mixed left-right 
handed sneutrino mediator. While this scenario was not constrained by previous data the Planck 
limits on N e g allow us to set a lower bound on the neutralino dark matter mass of 5.5 MeV. 



I. INTRODUCTION 

The mass and non-gravitational interactions of parti- 
cle dark matter remain unknown. Although many candi- 
dates have been proposed, particles that achieve the ob- 
served relic abundance through thermal freeze out have 
garnered the most attention [TH3]. While thermal relics 
are often assumed to have weak-scale masses this does 
not have to be the case. The thermal freeze-out mecha- 
nism works for a much broader range of masses and leaves 
room for models from keV to multi-TeV masses encom- 
passing the whole experimentally allowed range of cold 
dark matter. 

An upper bound can be set on the dark matter mass 
using partial wave-unitarity which leads to a mass limit of 
approximately 340 TeV (4| . Finding a robust lower bound 
is more difficult. General theoretical attempts have been 
made, for example by Hut and Lee- Weinberg [5J[B]. These 
can be avoided by keeping the dark matter particle in 
thermal equilibrium with, for instance, a light mediator, 
and a robust lower bound requires using experimental 
constraints. For example, a model independent lower 
bound for fermionic dark matter of m > 1 keV can be 
set by considering the phase-space distribution of dark 
matter in dwarf spheroidal galaxies [7] . Structure forma- 
tion indicates that dark matter cannot be hot so a similar 
bound of O(keV) can be placed by requiring that the free- 
streaming length is small [8]. A more model dependent 
lower bound, m > 10 GeV, can be placed by considering 
distortions due to energy injection by s-wave annihila- 
tion in the cosmic microwave background (CMB) [9rll4|. 
However this can be evaded if the annihilation is a p-wave 
process or if the dark matter annihilates to neutrinos. 
Bounds can also be placed from Big-Bang Nucleosynthe- 
sis |15[I16|. but these suffer from systematic uncertainties 
in the determination of the helium fraction Y p 17, 18 . 
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Here we find a new lower bound for cold thermal dark 
matter using data from the Planck satellite. The only 
light states in the Standard Model that it is possible to 
be in thermal equilibrium with are the neutrinos or elec- 
trons and photons, and so it is natural to assume that the 
dark matter will be coupled to these states. If this is the 
case and the dark matter decouples while non-relativistic, 
there will be a change in the effective number of neutrinos 
(N c g) through its effect on the neutrino to photon tem- 
perature ratio [i"9rf2T] . Using the results from Planck [22] , 
which has measured the CMB angular power spectrum 
with an unprecedented accuracy, we show that a bound 
of about an MeV can be placed on the dark matter mass, 
if the dark matter is in thermal equilibrium with the neu- 
trinos or the electrons and photons after neutrino decou- 
pling. The exact bound depends on whether the dark 
matter is a scalar or fermion and to which SM species it 
couples. The advantage of our method is that works even 
for a p-wave annihilation cross-section and annihilation 
to neutrinos - cases when limits such as energy injection 
into the CMB do not apply. Previous studies have used 
iV e fT to constrain the mass of dark matter annihilating 
exclusively into electrons |19j . Our analysis using the 
Planck data supersedes that study and we also consider 
the case of annihilation into neutrinos. 

These bounds apply to models such as [23TI26] and elec- 
tric dipole and anapole dark matter where the dark mat- 
ter annihilates into electrons or photons, and could be 
relevant to models which give a signal at electron re- 
coil direct detection experiments [2"Tl |2"8] . Fewer mod- 
els have been constructed where dark matter annihilates 
into neutrinos. As an application of the bound, we show 
how this can constrain a supersymmetric (SUSY) dark 
matter model where the dark matter is a light bino-like 
neutralino and mixed left-right handed sneutrino plays 
the role of the light mediator. Annihilation in this model 
is p-wave and the neutralino remains in equilibrium with 
the neutrinos. We show that a new lower bound can be 
placed on the neutralino mass of 5.5 MeV. 

In Section [H] we show how light dark matter annihi- 
lating into neutrinos or electrons and photons changes 
N e g and derive a bound on the masses of dark matter 
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coupling to electrons and to neutrinos. Following that in 



Sections III we construct a supersymmetric model with 
a dark matter neutralino and a light sneutrino which 
mediates neutralino annihilation into neutrinos. In the 
appendix we have collected together constraints on this 
model from invisible widths of the Z- and Higgs boson, 
rare meson decays, collider direct production and beam 
dump searches as well as constraints from SN1987A, large 
scale structure and direct detection experiments. 



II. A LOWER MASS BOUND FROM N cB 

In this section we show how a lower bound can be 
placed on the mass of cold thermal dark matter by using 
Planck's constraint on N Q g. This bound applies if the 
dark matter remains in thermal equilibrium with either 
photons and electrons or neutrinos until after the time 
when neutrinos decouple from electrons. 

The effective number of neutrinos, N c g, is a convenient 
way to parameterise the energy density apart from the 
photon contribution. It allows us to write the total en- 
ergy density as 



P = P~f 
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The factor (4/11) 4 / 3 is the fourth power of the neutrino- 
photon temperature ratio T®/T 7 = (4/11) 1 / 3 , assuming 
instantaneous decoupling of three neutrino species. In 
the standard cosmological model N e g — 3.046, reflecting 
that the actual neutrino-photon temperature relation is 
slightly higher due to decoupling effects GU] • 

The standard lore is that a value of N e g > 3.046 im- 
plies that there must be a new light species which is 
relativistic when the cosmic microwave photons decou- 
ple (at T ~ 1 eV). However, this is not the case and 
thermal dark matter with mass up to ~ 10 MeV can 
also change N e g, as we now argue using the example 
of electrons. In the standard cosmological model, the 
electrons do not directly contribute to the energy den- 
sity at recombination because they are non-relativistic 
at this time. However they do have an effect on A c ff 
because the electrons caused the photons to develop a 
higher temperature than the neutrinos. The definition 
of N e ff assumes that a definite temperature ratio holds, 
namely that T°/T 7 = (4/11) 1 / 3 . Had the electrons not 
been present, the photons and neutrinos would have had 
the same temperature and N e g would have been larger 
than three by a factor (11/4) 4 / 3 . 

In a similar fashion, thermal dark matter in ther- 
mal equilibrium with photons, electrons or neutrinos can 
change N e g by altering the neutrino-photon temperature 
ratio. In order for this to happen, the dark matter must 
transfer its entropy to either the electrons and photons 
or to the neutrinos, after the neutrinos decouple from 
the electrons at Tb ~ 2.3 MeV [3T] (we assume that new 



interactions do not change the neutrino decoupling tem- 
perature). If the dark matter is in thermal equilibrium 
with photons, electrons and neutrinos after neutrino de- 
coupling, it will reheat them equally and N c g will not 
change. Since the entropy transfer occurs when the parti- 
cle becomes non-relativistic, a change in temperature can 
happen if the mass m < few x Tb ~ 10 MeV. Cold ther- 
mal dark matter automatically satisfies the condition of 
being in thermal equilibrium until the temperature drops 
below its mass; typically, MeV-mass thermal dark matter 
drops out of chemical equilibrium when T ~ m/15. 

We now put the above argument on a more quantita- 
tive footing. We first consider the case where the dark 
matter is in thermal equilibrium with neutrinos. The 
case where it is in equilibrium with electrons then fol- 
lows in a straightforward fashion. The contribution of 
N v neutrinos and n particles (with mass m; and internal 
degrees of freedom gi) in thermal equilibrium with the 
neutrinos to the energy density is 



Pv+r, 

Pi 



i=l 



where 



I(x) = 



120 

7^4 
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ev±l 
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with limits I(oo) — and 1(0) = 1(8/7) for fermions 
(bosons) respectively. As usual, the sign + (— ) refers to 
fermion (boson) statistics. In general, T v ^ T®. There- 
fore, comparing eqs. JT]) and ([2|, we see that 
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Anticipating that the bound on mj is such that rrii 3> 
T v ~ 1 eV, we set I(mi/T v ) = so that 



(5) 



The ratio T v /Ty is determined by considering entropy 
conservation (see e.g. [20i 32 ). After neutrino decoupling 
at Td ~ 2.3 MeV, the entropy of the 'neutrino plasma' 
and 'electromagnetic plasma' are separately conserved so 
that (for T 7 < Td) 



T 



g* 



g*s:f 
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Here \ T indicates that g+ s should be evaluated at the 
neutrino decoupling temperature Td while and 
g* s:7 , defined through s„ = 2iT 2 g irS . M T^/A5 and s 7 = 
27r 2 g* s:7 T 3 /45 respectively, are the effective number of 
relativistic degrees of freedom in the neutrino and elec- 
tromagnetic plasmas. Explicitly, 
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where 
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with limits F(oo) = and F(0) = 1(8/7) for fermions 
(bosons) respectively and the sign + (— ) refers to fermion 
(boson) statistics. 

Again, anticipating that the bound on is such that 
rrii 3> T„ ~ 1 eV, we find that for particles only in ther- 
mal equilibrium with neutrinos, eq. ([5| simplifies to 



^-equil _ at— 1/3 
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For the case of particles in thermal equilibrium with 
electrons or photons, we again find eq. ^ and can use 
eq. ^ to find the new temperature ratio. In this case, 
we find 



7vr7,e-equil N 
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-4/3 



(10) 



where we have used F (to c /Td) ~ 1. 

Equations ([£]) and ( 10 ) are a central result of this paper 
as they allow us to put a bound on the dark matter mass 
by using the measured value of N e g from [22] 



N eS = 3.30IS si 95%; Planck+WP+highL+BAO, 

(11) 

which combines Planck precision measurement of the 
CMB, WMAP-9's polarisation data (WP) [33], SPT's 
high-£ measurement (highL) [34] and baryon acoustic os- 
cillations (BAO) measurements from large scale structure 
surveys [35H39"] . 

An extra (dark matter) particle in equilibrium with 
neutrinos thus increases iV e g. Requiring that N e g is con- 
sistent with eq. ( 11 ), we can exclude the following particle 



masses at 95% CL 



Real scalar 
Complex scalar 
Majorana fermion 
Dirac fermion 



to < 0.3 MeV 
to < 5.7 MeV 
to < 5.5 MeV 
to < 8.7 MeV 



(12) 
(13) 
(14) 
(15) 



Similarly, an extra particle in equilibrium with electrons 
decreases N c g. Thus, using the measurement eq. (Ill, 



we can exclude the following particle masses at 95% CL: 



Real scalar 
Complex scalar 
Majorana fermion 
Dirac fermion 



to < 4.2 MeV 
to < 7.6 MeV 
to < 7.5 MeV 
to < 10.3 MeV 



(16) 
(17) 

(18) 
(19) 



These bounds are independent of whether the annihila- 
tion cross section is s- or p-wave, but there are ways to 
evade them. For example, while dark matter annihilating 



into electrons would by itself decrease N e g, the presence 
of some extra relativistic degrees of freedom ('dark radi- 
ation') or sterile neutrinos would increase N e s [2T ] 140 ] . 
or the dark matter's abundance might arise from a non- 
thermal mechanism so that agreement with the Planck 
limit might be found. We do not consider such a scenar- 
ios in this article. 



III. APPLICATION: LIGHT NEUTRALINO 
DARK MATTER 



We now consider the application of the above bounds 
to a specific cold thermal dark matter candidate. In our 
model, the dark matter candidate is a light neutralino 
which remains in thermal equilibrium with neutrinos af- 
ter they decouple from the electrons. The annihilation 
cross-section is p-wave suppressed and the above bound 
provides the strongest constraint on the mass. 

The natural expectation for the mass of the lightest 
neutralino is about the weak scale, 0(100) GeV. However 
the absence of any signal of SUSY at the Large Hadron 
Collider (LHC) may indicate that the much-studied usual 
mechanisms of SUSY breaking are not realised, and/or 
that some tuning may be required. Yet SUSY with im- 
parity still has many desirable features such as gauge 
coupling unification and may still tame most of the clec- 
troweak hierarchy problem, and provides a flexible frame- 
work that is extremely useful to explore possibilities for 
the dark matter parameter space. In light of these ar- 
guments, we consider an extension of the MSSM with a 
spectrum which is dramatically different from those con- 
sidered in standard SUSY model building (summarised 
in fig.§. 

In the MSSM, the lightest neutralino is a superposition 
of bino (B), neutral wino (W°) and neutral higgsinos (h d 
and/jO): ' 



X° = N n B + N 12 W° + N 13 h° d + N 14 h° u , 



(20) 



where Nu are the components of the neutralino mix- 
ing matrix [41) . Searches at LEP impose a lower 
bound on M2, which, when the GUT relation Mi = 
5 tan 2 6*vyM2/3 w M2/2 is imposed, translate into a lower 
bound on M\\ in this case a lower bound of 46 GeV can 
be placed on the mass of the lightest neutralino jT5J [22] • 
However, when this assumption is relaxed and M\ and 
M% are treated as independent parameters, a massless 
neutralino is not ruled out by laboratory, astrophysical 
or cosmological limits if it is bino-like [33T|4"9"] . Therefore, 
here we consider the case when the dark matter is, for all 
purposes, purely bino-like. 

In this paper we remain agnostic as to the exact mech- 
anism that gives rise to a light neutralino. However, 
we briefly mention two mechanisms by which this could 
arise. First, a light bino-like neutralino can be obtained 
if there is a hierarchy between M\ and (M 2 , /i) , which can 
arise in hybrid models of supersymmetry breaking [50] . 
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In this case, in the limit M2,/i 3> mx,Mi, at tree-level 
the lightest neutralino mass is given by 



Mi - sin(2£) sin 2 V 



and the bino-component is 



liV 



11 



1 + sin 2 6w 1 " Z 



1-1/2 



1 - 0.01 



300 GcV 



The 
\Niz\ 



higgsino 
i sin^^T 



components 
and 



are 

and |JV U | w cos/V 1 ~ 
respectively, while the wino-component is negligible: 
N X 2^0{N 11 M 1 /M 2 ). 

Second, we can invoke a symmetry argument. In SUSY 
models with an i?-symmetry, gaugino soft-masses are for- 
bidden (Dirac masses are required for the wino and gluino 
but the bino does not require them). However, it is ex- 
pected that supergravity effects violate the (global) R- 
symmetry giving rise to a light bino- like neutralino, such 
as in |51j . In that article the Majorana neutralino mass 
is 



(21) 

(22) 
(23) 



sub-dominant: 

cos /9yl~ 



Mi 



11a 



47T COS 2 6 



9 x 10" 



m 3 / 2 



(24) 



and the neutralino is totally bino-like. If the gravitino 
mass is low, 0(100) GeV or less, then the lightest neu- 
tralino will be at the MeV scale. 



A. Thermal neutralino dark matter 

The arguments of Hut 5 and Lee- Weinberg [B] suggest 
that a neutralino lighter than 0(10) GeV cannot have the 
correct thermal relic abundance in the MSSM as all me- 
diators with Standard Model particles are 0(100) GeV. 
Indeed, detailed studies have borne this out [521454] . We 
will therefore need to consider an extension of the MSSM 
which has a new light mediator |55j in order to keep 
the neutralino in equilibrium with the thermal bath for 
longer. 1 The light neutralino can only annihilate to those 
states present in the early universe, namely electrons, 
photons and neutrinos. We choose as our mediator a 
light mixed left- and right-handed sneutrino, which keeps 
the neutralino in equilibrium with neutrinos through the 
diagrams shown in fig. [I] The mixing with a sterile right- 
handed component is necessary because the left-handed 
sneutrino couples strongly with the Z-boson. Constraints 
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FIG. 1: Feynman diagrams showing the interaction that keep 
the neutralino in chemical equilibrium with the neutrinos in 
the early universe. The mediator v\ is a light mixed left- and 
right-handed sneutrino. 



from the Z-width prevent it from being lighter than 
mz /2 (this is discussed in greater detail in the appendix) . 

Therefore, in addition to the MSSM matter content 
we introduce three right-handed neutrino superfields iVj , 
which are gauge singlets with respect to the Standard 
Model gauge symmetries. The additional super-potential 
term which arises is 



AW 



Y v H 



LiNn 



(25) 



where Li is left-handed leptonic chiral superfield, H u is 
the up-type Higgs superfield and for Sf7(2)-doublets, our 
notation is A ■ B = A 2 B 1 — A Y B 2 . When the scalar com- 
ponent of H u gets a vev (h u = vsin/3), where v = 174 
GeV, this gives rise to a Dirac mass for the neutrinos. 
Therefore, the magnitude of the entries Y are approxi- 
mately Yfi ~ 10~ 12 , and this matrix is diagonalised by 



the PMNS matrix [56H58] . We do not consider a Majo- 
rana mass term in this article. 2 

With the additional soft terms, the left- and right- 
handed sneutrino contribution to the scalar potential is 



Vaoft 3 rri~ L \v Ll \ £ + ml\hA 2 + A. t] h u ■ Li?ij + h.c. . (26) 



For simplicity we have assumed that m| and m| are 
generation independent and that Aij has real entries with 
the same flavour structure as Y^. Therefore, Aij is also 
diagonalised by the PMNS matrix. This means that af- 
ter diagonalising A^, the sneutrino eigenstates and 
n a have the same flavour structure as the neutrinos. For 



Of course, this not the only option for a light neutralino to obtain 
the observed relic abundance: it could occur via a non-thermal 
mechanism or some source of entropy could dilute the neutralino 
overabundance. 



2 A Majorana neutrino mass term could be generated radiatively 
along the lines of I59H65I by including a supersymmetric mass 
term m^NN and the accompanying soft B-term in the scalar 
potential. Wc do not investigate this further. 
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instance, the neutrino mass eigenstate Vi has approxi- 
mately equal electron, muon and tau flavour so v\ and 
n 2 also have approximately equal electron, muon and tau 
flavour. 

The diagonal A a terms induce mixing between the 
left- and right-handed eigenstates (without changing the 
flavour structure) so that the lightest mass eigenstate is 



sin 9\ i>1 + cos 9\ n° 



with mass 



to? cos 2 6>i 



9 • 2 

to~^ sin 



cos 26>i 



The mixing angle is given by 



tan 29; 



2AiV sin /? 



(27) 



(28) 



(29) 



TeV -- 



GeV -- 



MeV -- 



%2 



For definiteness, we assume that the flavour index is 
a = 2 for the lightest mass eigenstate so that it cou- 
ples to the electron, muon and tau neutrinos with equal 
strength (to a good approximation). In the remainder of 
the paper, we drop the subscript T' on the mixing angle 
of the lightest mass eigenstate. 

The left-handed sneutrino mass-term, tnp L = 
to| + | cos(2/3)to|, shares the soft-mass rh L with 
the left-handed selectron so that the ATLAS bound 
m SL > 195 GeV [Mj implies that m 9L > 200 GeV. 
Therefore, to have a light sneutrino requires the acci- 
dental relation to s w tan^m^. We assume that this 
cancellation occurs for the lightest sneutrino while the 
others remain at 0(100) GeV. To achieve such a mass 
requrires a tuning comparable to that required in the 
electroweak sector of the CMSSM [571155] . 

To summarise, the neutralino and sneutrino mass spec- 
trum is shown in fig. [2j In addition to the light bino- 
like neutralino, we also have a light mixed left- and 
right-handed sneutrino. We assume that all other sparti- 
cles, including the heavier sneutrinos and neutralinos are 
heavy enough to be consistent with direct search bounds 
from the LHC and LEP. 

We turn now to the calculation of the relic abundance. 
The dominant interaction contributing to neutralino an- 
nihilation with the neutrinos is the t- and u-channel ex- 
change of the mixed sneutrino shown in fig. [I] In the 
non-relativistic limit, t « u w — to 2 and the resulting 

p-wave thermally averaged annihilation cross section is 



7\V 



lA^n^sin 4 ^ 
(m\ +m| ; ) s 



x 



(30) 



103 pb {sin 9 

~b~T 



5 MeV 



35 MeVY 



(31) 



where as usual, x = to^o /T. Using the well known ana- 
lytic approximation for p-wave annihilation, for instance 



FIG. 2: A schematic diagram showing the structure of the 
neutralino and sneutrino sectors that we consider in this pa- 
per. As well as a sub-GeV bino-like neutralino, we have a light 
mixed left- and right-handed sneutrino. The other neutralinos 
and mixed sneutrinos have soft masses around the electroweak 
scale 0(100) GeV. We assume that the other sparticles are 
heavy enough to have evaded direct search constraints from 
the LHC and LEP. 



see [32] for details, we find that the thermal relic abun- 
dance is 



o oh 2 = . 15 ( * ) VSi 0*/15) 

Xl V2889.23 cm- 3 J g* s \ xt (cr /100pb) 



where do is defined through (a ann v) — aQX 1 and 



Xi = In 



0.076^5* 
I In (in 



1/2 



TOpim^ocio 
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0. 076gg* to p ito^o(To 



(32) 



(33) 



We have explicitly left in the dependence of today's en- 
tropy density 



s = 2889.23 cm" 



/_9*s_\ 
V 3.9091/ 



(34) 



because the light neutralino and light sneutrino lead to a 
change in the present day value of g+ s = 2 + 2\T v j (4T 7 ), 
through the change in the neutrino-photon temperature 
ratio. 

Figure [3] shows the sneutrino mass required to 
achieve the relic abundance il^ah 2 = 0.11 as a func- 
tion of the neutralino mass m^o for three different val- 
ues of sin 9. For sin 9 > 0.05, the required sneutrino 
mass is ~ few x to-^o . For sufficiently small val- 
ues of sin0, the required sneutrino mass is smaller than 
the neutralino mass so that the neutralino can no longer 
be the dark matter. For instance, for sin 9 = 0.04, the 
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FIG. 3: The mixed sneutrino mass ma 1 required to achieve 
ft^oh 2 = 0.11 as a function of the neutralino mass m ; o for 
three different values of sine. The upper dashed blue line 
is for sin6 = 0.1, the middle solid red line for sine = 0.06, 
and the bottom dash-dotted black line for sine = 0.04. The 



required sneutrino mass is m fl 



few x m^o. 

x i 



sneutrino mass (required to achieve the observed relic 
abundance) is approximately the same as the neutralino 
mass at ~ 19 MeV. In this article we do not consider the 
case where the mixed sneutrino is lighter than the neu- 
tralino (and thus, potentially, the dark matter particle), 
although we note that similar bounds deriving from N e g 
will also apply in that case. 



B. Constraining the model 

In the appendix we have collected together constraints 
on this model from invisible widths of the Z- and 
Higgs boson, rare meson decays, collider direct produc- 
tion and beam dump searches as well as constraints 
from SN1987A, large scale structure and direct detec- 
tion experiments. These constraints imply that the neu- 
tralino must be bino-like and the sneutrino mixing angle 
sin (9 < 0.11. However, these do not provide us with any 
constraint on the mass of the neutralino dark matter, 
apart from the O(keV) mass bound from free streaming 
and phase space arguments mentioned in the introduc- 
tion. 

Therefore, we now use the N e g bound from section |TT] 
to show that this provides the strongest mass constraint 
for neutralino dark matter in this model, improving the 
previous lower bound by three orders of magnitude. Since 
both the neutralino and sneutrino remain in thermal 
equilibrium with the neutrinos after they decouple, we 
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FIG. 4: 



The three lines show N e g as a function of the neu- 
tralino mass m x for three choices of sin 8. The sneutrino mass 
rrivx nas been set to the value in fig. plso that ^^oh 2 = 0.11. 
The green shaded area shows the 95% confidence region 
N e ff = 3.301q'5* from the Planck satellite in combination 
with WMAP-9's polarisation data, SPT's high-^ data and 
BAO measurements. We can exclude m^o > 5.5 MeV. The 
difference in the behaviour of the lines is due to the sneutrino, 
which contributes when m„ 1 < 10 MeV. 



can use eq. (|9| to show 



N cB = 3.046 



F(m x; /T D ) +F(m* 1 /r D )' 



3.046 



4/3 



In fig. |4jwe plot iV e g as a function of m^o for the same 
values orsm 9 as in fig. [3] The same colour coding is also 
used in this figure. The sneutrino mass has been 
fixed to the value shown in fig. [3] for a given value of 
m^o and sin 9. The shaded green band in fig. |4j shows 
the value of N c g measured by Planck's precision mea- 
surement of the CMB, together with WMAP-9's polari- 
sation data, SPT's high-£ data and BAO measurements 



at 95% CL, quoted in eq. (11). We find a lower bound 



> 5.5 MeV, consistent with the value quoted in 



eq. ( 14 1 for the case of M ajorana fermion dark matter. 
When the the sneutrino mass drops below ~ 10 MeV, it 
can also contribute to N c g. This explains why the three 
lines behave slightly differently and why the bound when 
sin 9 = 0.04 is slightly weaker. 

In order to use the constraint on N e g, we have 
implicitly assumed that the neutrinos decouple at 
Td = 2.3 MeV and the neutralino and sneutrino are only 
in thermal equilibrium with neutrinos. We now show 
that this is indeed the case. 

The neutralino and sneutrino do not lead to a new in- 
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teraction between electrons and neutrinos that is stronger 
than the weak interaction. Therefore, the only change in 
the decoupling temperature occurs because the expan- 
sion rate is slighter faster due to the presence of addi- 
tional relativistic particles [ID] . Neutrinos decouple when 
r ~ G^Tfj = l.Q&gi^T^/mpy. The extra neutralino and 
sneutrino slightly change Td because they contribute to 
g*. Therefore, the ratio of old and new decoupling tem- 
peratures equals the ratio of old and new values of gV % ■ 
This change is always very small. For instance, taking 
to^o = 5 MeV and = 20 MeV, we find the neutrino 
decoupling temperature increases to 1.02 Td- 

We now show that the neutralino and sneutrino are 
only in thermal equilibrium with neutrinos. In order for 
this to happen, we have to show that the neutralino and 
sneutrino are no longer in kinetic equilibrium with the 
electrons at 2.3 MeV. We first consider the neutralino 
elastically scattering with electrons. The rate of elastic 
scattering is 



Tel = (^x°e V ) n e 



(36) 



and the temperature at which the last elastic scattering 
is found from solving r o i = H. However, if the neu- 
tralino is non-relativistic, it drops out of kinetic equi- 
librium before this time because the typical momentum 
transfer in a single collision, Ap^o /p^o f 
small compared to the momentum of the neutralino 
Therefore, N ~ a/2/3 m%o /T collisions are required for 
to maintain kinetic equilibrium so, when the neutralino 
is non-relativistic, the temperature of kinetic decoupling 
is found from solving = H , where 



3/2T/m^o, is 



r 



kd 



3 T 



2 m^o 



(a^o e v)n e 



(37) 



In the limit that the electron mass is zero and left- 
and right-handed selectrons are degenerate the thermal 
averaged elastic scattering cross-section is 



(a^o e v) 



Gls^e w {El), ,:is; 

zi \ m, 



where (-E 2 ) — 12.9 T 2 . The neutralino drops out of ki- 
netic equilibrium with the electrons when 



T, 



kd-e 



2.3 MeV 



(lOMev) (93Gev) ' (39) 



With a more accurate numerical calculation that in- 
cludes a non-zero electron mass and a thermal average 
over Fermi-Dirac distributions for the electron and neu- 
tralino, we find that mg > 70 GeV for m^o = 10 MeV 
in order that T kd -c > 2.3 MeV. The XtLAS limit 
rrig > 195 GeV 66 implies that the neutralino drops out 
of kinetic equilibrium with the electrons well before the 
neutrinos decouple. 

The sneutrino is kept in equilibrium with the electrons 
through Z- and chargino-exchange. The chargino con- 
tribution can always be made small by increasing the 



chargino mass, so we only consider the Z-contribution. 
Ignoring the electron mass we find 



sin 4 OGl 



Ssin6» 4 l ,-4sin6l^ + l) (E 2 e 



(40) 

In this case, the sneutrino drops out of kinetic equilib- 
rium with the electrons when 



T kd _ c « 2.3 MeV 



(J^\ X, *(™L\. (41) 
V 20 MeV/ \sm6 J V ' 



Performing a numerical calculation with a non-zero elec- 
tron mass and a thermal average over Fermi-Dirac dis- 
tributions for the electron and sneutrino, we also find 



that sin 6» < 0.37 for 



= 20 MeV. This value of sin 6 



is well above the values we consider, which implies that 
the sneutrino, like the neutralino, drops out of kinetic 
equilibrium with the electrons well before the neutrinos 
decouple. 



IV. CONCLUSIONS 

In this article we have used the recent Planck measure- 
ment of N e g to set a lower bound on the mass of cold 
thermal dark matter. These bounds are summarised in 
eqs. [12] to [15] for the case when dark matter is kept in 
thermal equilibrium with neutrinos and eqs. 



16 to 19 for 



the case when dark matter is kept in thermal equilib- 
rium with electrons and photons. These bounds apply 
for both s- or p-wave annihilation cross-sections. These 
bounds apply to models already in the literature such 
as [23H26] where dark matter annihilates into electrons 
and could be relevant for models which give a signal at 
electron recoil direct detection experiments. 

The case where the dark matter remains in equilib- 
rium has received less attention and therefore we have 
constructed a supersymmetric model to realise this sce- 
nario. In this model the dark matter is a light bino-like 
neutralino and it is kept in equilibrium with neutrinos via 
a mixed left-right handed sneutrino mediator (see fig. [2| . 
In the presence of a such a mediator it is possible to 
achieve the correct relic density (see fig. [3]), leading to 
changes in N c g (see fig. [4]). The Planck measurement of 
N e g leads to a lower bound of 5.5 MeV on the neutralino 
mass in this scenario. 

If the dark matter couples to electrons and there is also 
some "dark radiation" such as sterile neutrinos present 
the bounds may be weakened. Our bounds can also be 
evaded if the dark matter's abundance has a non-thermal 
origin (for example asymmetric dark matter), or if the 
dark matter is in thermal equilibrium with electrons, neu- 
trinos and photons in which case there is no change in 
the standard neutrino to photon temperature ratio. 
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Group. The neutralino contribution to the decay width 
is 



16tt 





a/ 4m| ; \ 








V m l ) 



3/2 



(A.5) 



Appendix: Constraints on a light neutralino and a 
light mixed sneutrino 

In this section we will show that the scenario we 
have outlined satisfies current constraints from particle 
physics experiments and astrophysical and cosmological 
processes. The light bino case has already been consid- 
ered extensively in |47H49] where it was shown that cur- 
rent constraints are satisfied. We update this analysis 
to include invisible Higgs decays to neutralinos, and per- 
form a qualitatively similar analysis for the light sneu- 
trino case. We find that all bounds can be evaded. The 
most stringent bound comes from the invisible Higgs- 
width which puts a limit of sin# < 0.11. 



Z- and h,-boson widths 



withy' 1 *?*? = g(N 12 - tan%iVn)(sino;iV 1 3 + cosaWi 4 ). 
In the decoupling limit and using N 13 /Ni 4 « — tan/3, 
valid when the neutralino is very light, we have, 



(jpm z sin 



%^/2-k 



m h sm 2 2/3N^(l-N^) 



79 MeV x sin 2 2/3^(1 - N^) 



(A.6) 
(A.7) 



For any value of Nu, tan/3 can be chosen so that the 
contribution to the higgs width is below the experimental 
limit. For instance, if Nu > 0.9(0.99), the invisible width 
is below the limit if tan /3 > 7.6(2.4). 

One mixed left- right-handed sneutrino complex scalar 
contributes 



Y 2 



1 - 



(A. 



The light neutralino and sneutrino couple to the Z- 
and higgs-bosons. We consider the contribution to the 
partial width of each in turn. Considering the Z first, 
the neutralino contribution is 



12V2-K 



'N 



166 MeV x cos 2 2/3 (l 



handed sneutrino complex scalar is 



sin 



G 



F m z 



24tt y/2 
83 MeV x sin 4 9 



1 - 



4m,-,j 



m r 



4m 2 - \ 


3/2 






m z J 






(A.I) 




(A.2) 


left-and 


right- 






^3/2 


(A.3) 




(A.4) 



< 3 MeV [70], we have that N n 



> 



With the limit T z ^inv 
0.87 and sin0 < 0.44, far from the parameter range we 
consider. 

The invisible width of the Higgs is a notoriously diffi- 
cult quantity to constrain at hadron colliders as the total 
width is quite inaccessible. Without some assumption 
about the Higgs couplings it is not possible to extract a 
meaningful constraint on this quantity. However, if we 
assume that the Higgs couplings are as given by the SM 
and fit to current Higgs data while including an invisible 
contribution to the wdith, the invisible branching ratio 
is constrained to be less than 15-20% [7lTf74"j or equiva- 
lently T h ^ inv < 0.82 MeV for an SM width of 4.1 MeV, 
as suggested by the LHC Higgs Cross-Section Working 



where 



Ai 
2V2 



cos a sin 28 



rriyy sin(a + j3) sin 2 

y/2v COS 2 W 



(A.9) 

The first term comes from Ah u i>n* -term while the sec- 
ond term is the usual contribution of the (left-handed) 
sneutrino. In the decoupling limit and using the relation 
Ai^sin/3 = tan 6*(m? — mf ) 



tan 0m^ L , we have 



Y, 



so that 



sin 2 BM\ 
V2v 



1 - 



0.82 MeV x 



mw 



sin 




cos 2/3 
cos 2 6 



w 



(A.10) 



1 + 0.21 



/ 195 GeV \ 
V m *L J 



195 GeV 



cos 2^ 

-0.98 



V 



1.21 



(All) 



The contribution to the higgs width is large because of 
the large A-term, required to have reasonably large left- 
right-mixing (c.f. eq. (291). However, we can always 
choose parameters such that we are below the current 
experimental limit. 



2. Meson decays 

We next consider constraints from partial widths of 
invisible meson decays. We briefly summarise the neu- 
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tralino case, previously investigated in before dis- 
cussing the implications of their results for the sneu- 
trino. The authors of [49] find that lower bounds cannot 
be set on the neutralino mass in minimal flavour viola- 
tion (MFV) scenarios. Only in the case of non-minimal 
flavour violation can the branching ratios be enhanced, 
a situation we do not consider in this article. Two-body 
decays of pseudoscalar mesons to neutralinos proceed at 
tree-level via an exchanged squark. Two-body decays 
to sneutrinos are loop-suppressed relative to this, and so 
pseudoscalar decays do not constrain the sneutrino mass. 

One can also consider three-body decays such as K~ — > 
7T — />i z>i, where the Standard Model process K~ — > -k~vi> 
are already suppressed. These have also been calculated 
for light neutralinos in [49] . where it was found that the 
ratio of the SUSY contribution to the observed branch- 
ing ratio was at most 10 -6 . While reproducing the full 
results of [55] is beyond the scope of the paper, we see 
no reason why similar results for the sneutrino final state 
should not also be true. 



3. Mono jet and monophoton searches 

An active area of recent research has been using col- 
lider data to set limits on theories of dark matter. Of 
particular interest to us are the studies [751 175] , the first 
of which deals with mono jet limits from the Tevatron 
and projected limits from LHC, and the second with 
monophoton searches from LEP. We consider each in 
turn. 

The authors of [75] adopt an effective field theory ap- 
proach, writing down higher dimensional operators which 
couple two dark matter fields to two quarks or gluons, 
suppressed by a scale M*. For light dark matter, the 
bounds they derive on are 0(10) GeV for Tevatron 
results, with projected limits of 0(100) GeV at LHC at 
14 TeV. As the main mediator in the MSSM for this pro- 
cess are squarks, direct search limits are already far above 
these bounds. The presence of a light sneutrino does not 
affect these limits, as processes like pp — ¥ v\V\ are loop 
suppressed relative to pp — > XiXi- 

At LEP there are two processes of interest leading to 
the monophoton final state. These have been considered 
in [75] in an effective field theory context, setting lim- 
its on A the suppression scale of the higher dimensional 
operator. For the process e + e~ — > XiXil this is approxi- 
mately A ~ mg cos Qw / V 47ra leading to > 120 GeV. 
The other process which contributes to the monopho- 
ton final state is e + e~ — > v\V\^ via a t-channel chargino 
(there is also a Z-mediated contribution but this is sub- 
dominate when the chargino is light). Here the sup- 
pression is approximately A ~ sin 6w / (V 47ra sin 9) : 
with sin# = 0.1 we find m^± > 21 GeV. These bounds 
are easily satisfied. 



4. Beam-dump searches 

It has been proposed that high-luminosity fixed-target 
neutrino experiments can be sensitive to sub-GeV dark 
matter [77 . In order to assess whether we are constrained 
by them we adapt the calculation of Ref. [75] to our case. 
They investigate sub-GeV dark matter models with both 
a vector and scalar portal, the scalar portal being the 
case relevant for us They find that the sensitivities of 
MINOS, T2K and MiniBooNE are too weak in the scalar 
case to set limits (unlike the vector case). In addition to 
the Higgs mediated diagram contributing to p\\ — > pxl 
scattering which they consider, we should also include 
the contributions from squark exchange. However, as 
one expects the squark mass to be ~ 0(1) TeV 
this contribution should be suppressed by approximately 
m t,/ m t ~ 10~ 4 . Our model is thus unconstrained by 
beam-dump searches. 



5. Bounds from SN1987A 

The neutrinos detected from SN1987A were in good 
agreement with the theoretical expectation. Requiring 
that new particles or new interactions do not significantly 
change the neutrino burst allows bounds to be placed 
on a wide variety of new particles and interactions (see 

e.g. eg). 

The scenario we consider here, namely, of MeV-mass 
particles with stronger-than-weak interactions with neu- 
trinos and very weak interactions with electrons and nu- 
cleons, has not been specifically considered in the litera- 
ture. Bounds have previously been set on new particles 
that have stronger-than-weak interactions with nucleons 
(or electrons) in addition to stronger-than-weak inter- 
actions with neutrinos, leading to a constraint ofm > 
10 MeV [80]. In our case, the interaction with electrons 
and nucleons is 'weaker-than-weak' so these bounds are 
not applicable. Our scenario shares more similarity with 
the Majoron [81H83J . a Goldstone boson that has large 
interactions with only the neutrinos, so we adapt bounds 
placed on that scenario. 

We first consider the energy loss argument (see e.g. [54T - 
186] for the Majoron case). The new particles can act as an 
efficient energy loss mechanism if they are so weakly in- 
teracting that they freely escape from the supernova core, 
resulting in a shorter Kelvin-Hclmholtz cooling phase. 
Here, we estimate the mean free path of the neutralino 
to determine if it escapes freely. We have numerically 
calculated the thermal average of the scattering cross 
section (assuming Fermi-Dirac distributions for the neu- 
trino and neutralino with zero chemical potential). For 
instance, for m x = 1 - 20 MeV, m 0l = 10 - 100 MeV 
and sin 6 = 0.1, we find (a v ^ov) rj 10~ 30 cm 2 . This 
cross-section is large because the sneutrino is on-shell. 
Assuming a neutrino density n v ss 10 33 cm -3 and a core 
temperature T = 10 MeV [57] we find that the mean free 
path of the neutralino is ~ 10~ 3 cm (this result is not 
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particularly sensitive to the choice of core temperature). 
A similar calculation holds for the sneutrino. Since the 
mean free path is so much smaller than the geometric 
dimension of the supernova (~ 10 km), the cooling from 
neutralinos and sneutrinos is very inefficient and should 
not decrease the Kelvin- Helmholtz cooling phase [79]. 3 
It has also been argued that the large interaction rate 
of the neutrinos with thermal neutralinos and sneutri- 
nos increases the cooling time because it takes longer for 
the neutrinos to diffuse from the core to the neutrino- 
sphere [58]. However, a more careful analysis treating 
the neutrinos as a relativistic Fermi gas indicate that 
this argument is false [59"] . 

A more robust bound can be placed considering the sig- 
nal dispersion of v e . The detection of rVs precludes the 
possibility that they were removed due to 'secret interac- 
tions' with dark matter particles as they journeyed from 
the Large Magellanic Cloud (LMC) to Earth [9Ql[9T]. In- 
tegrating the dark matter density along the line of sight 
between the LMC and the Earth and assuming the Milky 
Way's dark matter density follows p — po{lo/l) 2 with 
po = 0.3 GeV cm -3 and Iq — 8 kpc, we can put a bound 
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" 25 cm 2 



(A.12) 



on the dark matter scattering cross-section by requiring 
that that the mean free path for neutrinos is longer than 
the LMC-Earth distance. 

In the regime where > E v > m^o the elastic scat- 
tering cross-section is 



n -\N I 4 ™ 4 ft G 2 F m% S m i e w m^E, 



7/7 



37T 



10" 34 cm 2 



1 MeV/ 

E„ \ /sintA 4 /20 MeV x 4 



10 MeV 



0.1 



(A.13) 



(A.14) 



well below the limit in eq. (A.12) for the parameter range 
we consider. 



dark matter thermal decoupling. One can identify three 
distinct scenarios (93) . namely: scenario A defined by 
I\,_ e > T u _ dm and T dm - V <T V ~ r„_ e , scenario B de- 
fined by L_ e < Tv-dm and Ydm-u < IV ~ Tv-dm and 
finally scenario C, defined by Tdm-v > IV- We use 
to denote the interaction rate between species i and j, 
and Y v to denote the total interaction rate of the neutri- 
nos. 

Scenario B assumes a much larger dark matter- 
neutrino cross section than the neutrino-electron elas- 
tic scattering cross section at T ~ 1 MeV and implies 
that the neutrino thermal decoupling is actually modi- 
fied. The effect of such interactions has been studied in 
detail in [SU] H01 US] . Scenario C requires that the dark 
matter is coupled to free-streaming neutrinos and there- 
fore assumes very specific values of the dm-neutrino cross 
section, as discussed in [M] in the case of MeV dark mat- 
ter. Finally scenario A corresponds to very small cross 
sections (much smaller than the neutrino-electron cross 
section at 1 MeV) and the standard value of the temper- 
ature for the neutrino thermal decoupling. 

One means of parametrising this is by using the 
opacity Q = {<J xu v) /m x , the thermal average of the 
neutrino-dark matter elastic cross-section divided by 
the dark matter mass [SHI 155] . The Sloan Digital 
Sky Survey (SDSS) [H] constrains this to be Q < 
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-42 



r 2 MeV -1 



In the extreme non-relativistic limit 



where mo 1 > m^o 3> E v , we find 



Q 



3G 2 F m% sin 4 W | jV n | 4 sin 4 8(E 2 ) 
2tt m^o (to|o -mlj 2 

3 x 1(T 53 cm 2 MeV -1 

'sin(9\ 4 / 1 MeV' 



0.1 




(A.15) 



(A.16) 



where we used (E v ) « 12.9 T 2 with T v - 0.75 T lQ . Since 
this value corresponds to a very small cross section for 
MeV neutralinos, our model corresponds to the scenario 
A and there is no significant damping of the matter pri- 
mordial fluctuations. 



6. Bounds from large scale structure 

Dark matter interactions with Standard Model par- 
ticles can wash out primordial matter fluctuations and 
induce a damping of the matter power spectrum at 
large cosmological scales [92 94 . In the case of dark 
matter-neutrino interactions, the magnitude of the ef- 
fect depends on whether or not the neutrinos have en- 
tered their free-streaming regime at the moment of the 



3 Since the neutralinos do not freely escape the supernova, the 
SN1987A bounds on the selectron and squark masses from [47| 
do not apply. 



7. Direct detection 

As is well known, dark matter particles with sub- 
GeV mass do not have enough kinetic energy to cause 
a nucleus to recoil with an energy above the low-energy 
threshold of conventional dark matter direct detection 
experiments. However, in the sub-GeV mass range, dark 
matter particles can cause single-electron ionisation sig- 
nals by scattering off electrons and this signal is de- 
tectable [271 OH]- The first limit on the electron-dark 
matter scattering cross-section comes from 15 kg-days 
of data collected by the XENON10 experiment. They 
constrain the scattering cross-section to be less than 
a e = 6 x 10 -36 cm 2 for dark matter of mass 10 MeV. 
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Considering \i + e Xi + e an d assuming degenerate 
left- and right-handed selectrons, we find 

^ = S(S) 4G|sin4 ^ (AJ7) 

« 3 x 10- cm2 f^x y , (A . 18) 



where ^f e x? i s the neutralino-electron reduced mass. This 
is far below the XENON10 limit and below the projected 
limit er e ~ 1CP 43 cm 2 for a germanium detector with 
1 kg- year of data [57]. Therefore, it is unlikely that this 
model will be testable at future direct detection experi- 
ments. 
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